Charge and spin dynamics of interacting Fermions in a one dimensional harmonic trap 
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We study an atomic Fermi gas interacting through repulsive contact forces in a one dimensional 
harmonic trap. Bethe-Ansatz solutions lead to an inhomogeneous Tomonaga-Luttinger model for 
the low energy excitations. The equations of motion for charge and spin density waves are analyzed 
both near the trap center and near the trap edges. While the center shows conventional spin-charge 
separation the edges cause a giant increase of the separation between these modes. 

PACS numbers: 73.22.Lp, 32.80.Pj 
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Recent advances in cooling technology have allowed to 
reach the degenerate regime of Fermionic quantum gases. 
Different hyperfine states effectively correspond to spin 
polarizations of spin- ~ particles of respective densities n-f 
and nj, 0. Furthermore, the spatial dimensionality of the 
gases can be reduced using either a hollow beam setup Q 
or arrays of microtraps [3j so that one dimensional (ID) 
gases can be studied, where interaction effects are known 
to be most pronounced. As a consequence of pure s-wave 
scattering in 3D Q the interaction is short ranged || and 
acts only between particles of opposite spins. Also, the 
strength of the forces can, in principle, be varied over 
wide ranges |6| by tuning the Fcshbach resonance 0- An 
optical lattice along the trap can further enhance corre- 
lations by reducing the bandwidth and thus the kinetic 
energy H,@. 

As opposed to higher dimensions, total charge (p = 
n-f + ?i jj and relative spin (a — nj — nj,) density waves 
comprise the only low energy excitations in ID; no 
Fermionic quasiparticies can be excited in this reg ime. 
The Tomonaga-Luttinger liquid (TLL) theory 0, [n| 
describes the properties of homogeneous ID systems en- 
tailing that charge and spin modes move at different ve- 
locities when interactions are present. This spin-charge 
separation is considered as a hallmark of TLL behavior. 
A parabolic trap potential, however, causes the particle 
density to vary along the trap. To treat such an inho- 
mogeneous gas cloud, new Boson representations of the 
Fermi operators have been introduced Here, we fol- 
low the other route put forward recently by Recati et 
al. Q and consider an inhomogeneous TLL with IE- 
dependent parameters, assuming a trap potential which 
is slowly varying on the scale of the Fermi wavelength. 

In the bulk of the gas cloud, away from the trap edges, 
interactions are sufficiently weak to justify perturbative 
estimates to the (local) TLL model parameters. In this 
regime, the picture of charge and s pin waves moving at 
different velocities is recovered 0, Il4j . Near the trap 
edges, however, the gas density decreases and interac- 
tions become very strong. In this regime we employ the 



Bethe-Ansatz solution for spin-i particles with contact 
forces [l5L Hil Il7| to compute the TLL parameters ex- 
actly |18l [l9j . To leading order in the inverse interac- 
tion strength we can access the dynamics analytically. 
Charge modes are found to propagate in a similar way 
as in the bulk of the gas cloud. At the edges they are 
reflected and thus keep oscillating in the trap until damp- 
ing processes become significant. Remarkably, spin den- 
sity waves show an exponential slowing down of their ve- 
locity and, ultimately, accumulation at an edge without 
reflection. Thus, Fermionic ID quantum gases in shal- 
low confinements should exhibit a giant increase of spin- 
charge separation, much more pronounced than electrons 
in quantum wires |2fl | . The effect of spin accumulation 
should be detectable: one way would be to observe (e.g. 
by fluorescence measurements) the time evolution of an 
initial spin-up density peak containing charge and spin 
modes of equal amounts, pi = (p + a)/2. 

Without parabolic confinement the quantum gas is de- 
scribed by the Hamiltonian (h = 1) 
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Effectively, by virtue of the Pauli principle, only Fermions 
of opposite spins are interacting. As long as the 3D scat- 
tering length a is much smaller than the transversal width 
d of the gas cloud, the interaction strength c = 2a/d 2 

For given particle density n, the ground state energy 
density £hom(^) is obtained exactly from the Hamilto- 
nian via the Bethe-Ansatz [l^. In presence of the 
longitudinal trap potential V(x), we determine the in- 
homogeneous profile of the particle density n(x) in the 
spirit of the local density approximation (LDA) || by 



minimizing 



E[n] — J dx [eham(n(x)) + V(x)n(x) — p,n{x)] 



(2) 



This is adequate for slow variations of V compared to the 
Fermi wavelength and for large particle numbers, N 3> 1, 
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FIG. 1: LDA density profile of interacting Fermions in a trap 
of length 2R for dimensionless interaction strength u — 0.3, 
using the full Bethe-Ansatz equations numerically (solid line) , 



Eq. JSJ (dashed line), and Eq. (dotted line); 



no 



-moJT-R 



is the density of a non-interacting Fermi gas at the same fi. 



where we can ignore Friedel oscillations occurring on the 
scale of the Fermi wavelength Tr/2kp The chemical 
potential u ensures that N — J dx n(x) and determines 
the length 2R of the gas cloud through V(R) — fi, inde- 
pendent of V(\x\ < R) and the form of ehom(^), provided 
d n £houi{n = 0) = |21|. Assuming now a harmonic trap 



potential V(x) 
Eq. J2J 



lo^x of frequency wt, one finds from 



d n £hom(n(x)) + —lo t x - fi = 



(3) 



The resulting density profile n(x) is depicted in Fig. 

To proceed we now focus on two regions of particular 
interest. Near the center of the trap, where the density 
n is large enough to validate perturbative expressions in 
the interaction strength, one finds to leading order in c/n 

M 
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Inserting Q into © gives a modified Thomas-Fermi pro- 
file 



n pei t(0 = —vfo ( v 1 



u 2 — £ 2 — u 



(5) 



Here we have introduced the dimensionless coordinate 
£ = x/R and interaction strength u = c/mnvFo where 
Vfo = is the Fermi velocity of a homogeneous, non- 
interacting system at chemical potential fi. On the other 
hand, near the trap edges, c > n so that we can eval- 
uate the Bethe-Ansatz solution to leading order in n/c, 
yielding 4 8j 



n 2 n 3 ( An 

EedgeW = -g— II lli^ 



(6) 



Inserting JJJ into J2J gives 
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(7) 



As seen in Fig. ^ the analytic forms (J5J and {7J describe 
the trap density accurately in their regions of validity. 
We stress that n e dge(£) vanishes with infinite slope at the 
edges £ — > ±1, which is not seen from the perturbative 
expression n per t(0- 

The inhomogeneous TLL Hamiltonian is of the form 



H 



TLL 



v—p,a 



J2 / dx[vj„(x)(eUx)f+vMx)(K(x)) 2 



(8) 

where primes denote spatial derivatives, and 0„ and $„ 
are the usual Bosonic phase fields in the charge (y = p) 
and spin {y = a) sectors, obeying [<& v (x),®' v ,(x')\ — 
i5 uv i5{x — x'). They define density excitations v — 
(& v )/^/tt. The velocity parameters vmAx) follow from 
exact compressibility relations 0, H3, 13 



vt$ v = — 



2 9 2 ehom 



7T dv 2 



(9) 



while symmetries fix the values of velocity parameters 
vjp — 7ro/2m owing to Galilei invariance and vj a — 
«Nct as a consequence of SU(2) invariance in the spin 
sector. The model ||SJ| can be justified if v <§; n which, 
evidently, fails too close to the trap edges where n — > 0, 
see below. 

After separating off the time dependence ~ e lujt in the 
Heisenberg equations of motion, the eigenvalue equations 
for the density excitations become 



- VKuVlvv" - (2Vn„UJk + VNvV'jv) V 



(10) 
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As boundary conditions we require that charge and spin 
currents vanish at the trap edges, j v (x = ±i?) = 0. In 
view of (JBJ and the identity j„ = 0„/-y/7r, this boundary 
condition becomes vj v v' Nv v + v j v vhvv' = at £ = ±1 
which, together with (|10|l . governs the dynamics of den- 
sity wave packets. In the homogeneous case ^(n/j)i/ 
and v'qj ,js vanish and l|l(J|l simplifies to a wave equa- 
tion. In the inhomogeneous case, without interactions, 
all four TLL velocity parameters are equal and given by 
v f(0 = — £ 2 which follows from Eqs. © and (0 

for ehom(«) = 7r 2 n 3 /24m. Then the solutions of (fTTTf) 25] 



1 



■. cos (u) arccos £) 



(11) 



with integer to — w/wt, constitute a complete orthonor- 
mal set with regard to the measure — ^|J=5 [jfl |. 

Near the trap center, we obtain from Eqs. (0) and © 
to leading order in u 



V F = V F0 



(vT 



(12) 
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V-Np = v F + UV F0 



(13) by 



uv F0 ) 



v^y/l-x^/Rl , (14) 



where v^o = upoVl — 3u and = Ry/l — 3u/2. The 
solutions of I|1U[1 now become 



Pn{x) = 

a n (x) = 
•<(<*) 



n/2 



^(l + u 2 )R 2 -x 
1 



= cos n arccos — , (16) 

VW^x 1 V R*} 



where the C„ denote ultraspherical polynomials |27| . 
i?2 = flyl — u/2, ft = u/(2 — u), and the eigen- 
frequencies are cJ P)n = WT\/n(" + + ^)) w o-,« = 

wt«V (1 — 3u)(l — 3u/2). Thus, in the charge sector 
the spectrum is no more equidistant (see also Refs. 0, 
EH]), although the lowest mode (n — 1) remains unaf- 
fected by interactions in accordance with Kohn's theo- 
rem [2^. Eqs. (|15fl and l|16f) describe how weak interac- 
tions smoothly deform the non- interacting solution (|llfl ; 
note that c| 0) (£) = ^ cos( n arccos £) . The inequalities 
R2 UJ p,n > nvpQ and R a Lo an < nv F o reflect the enhanced 
and suppressed dynamics of charges and spins, respec- 
tively. 

Near the trap edges, in the non-perturbative regime 
where VR 2 ~ x 2 < uR (but still R - \x\ > fcp 1 ), inter- 
actions alter the above picture, even qualitatively. Using 
Eqs. jnj and 0, and employing the Bethe-Ansatz re- 
sult by Coll [lflj . we get from to leading order in 
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(17) 
(18) 

(19) 



where /3 = 8 In 2/ (Sir 2 u). Disregarding for the moment 
subleading terms oc /3 in Eqs. l(T7|l and ifTHl) . we see that 
the charge sector exponent K p = vj p /v^ p approaches 
the value 1/2 close to the edges, which is consistent with 
the limiting value of K p in the Hubbard model when the 
filling goes to zero [29j. The plasmon velocity y/vWpVjp, 
on the other hand, takes exactly the value of the non- 
interacting system [cf. Eqs. I|12fl and (|13l) for u — ► 0]. 
This amazing property originates from the fact that the 
particle density is reduced by a factor 1/2 at infinite in- 
teractions, compared to the non-interacting system. Fur- 
thermore, since Eq. IjlOl) only contains products vnuVjl, 
and their spatial derivatives, we conclude that even the 
whole charge dynamics coincides exactly with the non- 
interacting dynamics l|ll|) sufficiently close to the trap 
edges. 

To leading non-trivial order in /3, Eq. (|1(J|) can be trans- 
formed into a damped Mathieu equation, that is solved 



p n (x) 



VR 2 - x 2 



se2n(q, \ arcsin£) , nodd 1 
ce2«((?, \ arcsin£) , neven J 
(20) 

with the Mathieu functions se n and ce„. The eigenen- 
ergies u> n follow from the characteristic values [27| 
ainisLn) = ^lu 2 — 9/3 2 of the Mathieu functions while 
q n = -f3{2Q 2 n - 1). When 0^/1 - £ 2 -> 0, the p n (x) turn 
into the non-interacting solutions i|ll|l , in agreement with 
the above observations for the vicinity of the trap edges. 
In the spin sector the corresponding solutions, 



cr(x, lo) 



R 



R 2 



sin 1 / 3utu , , . , 
, , \ Atanh^ , 21 

r ^ 1 COS J \ Wt 



differ qualitatively from the non-interacting solutions 
(| 1 1 P . as well as from the weakly interacting ones Hl(j|l . 
They are rapidly oscillating near the trap edges which 
manifests extremely slow spatial propagation of spin den- 
sity waves, as we discuss now. 

With the eigenfunctions for charge and spin densities 
at hand, we are in the position to solve for the time evo- 
lution of an initial wave packet vq{x) in the trap center, 
or near the edge. After preparation, the packet splits into 
left- and right moving parts, v + and and the charge 
and spin constituents separate. Near the edge the time 
evolution follows 



p±{x,t) 
x' 

a±(x,t) 



R 2 



R 2 



-Po(x') 



R cos (arccos £ ± us^t) 



1R 2 



-a Q (x') 



2 R 2 

R tanh ( Atanh £ ± — t 
\ 3m 



(22) 



(23) 



Subleading corrections of order (3 require a summation 
over the eigenfunctions (|20|l and can no longer be ex- 
pressed in the simple form (|22|) . 

It follows from Eq. i|22|) that the charge density ex- 
hibits a temporal evolution that roughly resembles that 
in the absence of interactions: the right moving part of 
the initial wave packet slows down somewhat due to the 
decreasing Fermi velocity, and then is reflected at the 
edge (cf. Fig. |2^)- We thus expect charge density exci- 
tations to keep oscillating in the trap, before damping 
mechanisms set in. On the contrary, the right moving 
part of the spin density wave is slowed down even ex- 
ponentially fast, see Eq. I|23l) . and is not reflected (cf. 
Fig. |2h) ■ We mention that for long times Eq. gives 
<j±(x, t — > 00) — > A 5(x ± R) where A = J dx cr (x) is 
the conserved magnitude of the initial spin wave packet. 
This follows from the finite slope of v a at £ = ±1, as op- 
posed to the infinite slope Up( N /j) at £ = ±1 in the charge 
sector. However, our analysis is only valid up to times 
of order t r w — In 



2m Ru T /i _ c2 



3 In jtj- 
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both in the bulk and near the edge of the trap. 

We thank R. Egger, A. Komnik and W. Zwerger for 
valuable comments. 



FIG. 2: (Color online) Propagation of a charge density (a) 
and a spin density (b) excitation for dimensionless interaction 
strength u — 0.3. Time steps are AtuiT = 0.25 as indicated. 
The charge wave is reflected at the edge of the trap, similar 
to the non-interacting case, while the spin wave exponentially 
slows down without reflection. 



then a±(x,t > t c ) exceeds n(x). This time is large for 
small initial displacements <To so that a dramatic separa- 
tion of the spin and charge peaks will have occured well 
before effects on the length scale of the Fermi wavelength, 
that have been disregarded here, become relevant. 

In conclusion, we have studied interacting Fermionic 
atoms of two spin species in an effectively one dimen- 
sional harmonic trap by exploiting solvability by the 
Bethe-Ansatz for contact forces. We have evaluated the 
particle density profile within the local density approxi- 
mation. Near the trap center we confirm the occurrence 
of spin-charge separation. Near the edges of the trap in- 
teractions affect the dynamics in an unexpectedly drastic 
way. While in the immediate vicinity of an edge charge 
density waves move as if they were non-interacting with 
reflection at the edge, spin density waves are not reflected 
and accumulate at the edge. This establishes an even 
more pronounced spin-charge separation than in homo- 
geneous systems. Experimentally one could confirm our 
predictions e.g. by selectively evaporating parts of the 
'up' spin component p-f in a small region of space and 
observe the right and left moving spin and charge density 
waves separating from one another (since 2p| = p + a), 
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